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Abstract We niodify the Einstein-Schrodinger theory to include a cosmolog- 
ical constant which multiplies the symmetric metric, and we show how the 
theory can be easily coupled to additional fields. The cosmological constant 
Az is assumed to be nearly cancelled by Schrodinger's cosmological constant 
All which multiplies the nonsymmetric fundamental tensor, such that the to- 
tal A = A^+Ai, matches measurement. The resulting theory becomes exactly 
Einstein-Maxwell theory in the limit as ^oo. For \Az \ ~ 1/ (Planck length)^ 
the field equations match the ordinary Einstein and Maxwell equations ex- 
cept for extra terms which are < 10~^^ of the usual terms for worst-case field 
strengths and rates-of-change accessible to measurement. Additional fields can 
be included in the Lagrangian, and these fields may couple to the symmetric 
metric and the electromagnetic vector potential, just as in Einstein-Maxwell 
theory. The ordinary Lorentz force equation is obtained by taking the di- 
vergence of the Einstein equations when sources are included. The Einstein- 
Infeld-Hoffmann (EIH) equations of motion match the equations of motion for 
Einstein-Maxwell theory to Newtonian/Coulombian order, which proves the 
existence of a Lorentz force without requiring sources. This fixes a problem 
of the original Einstein-Schrodinger theory, which failed to predict a Lorentz 
force. An exact charged solution matches the Reissner-Nordstrom solution ex- 
cept for additional terms which are ^ 10^^^ of the usual terms for worst-case 
radii accessible to measurement. An exact electromagnetic plane- wave solution 
is identical to its counterpart in Einstein-Maxwell theory. 
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1 Introduction 

The Einstcin-Schrodinger theory is a generaUzation of vacuum general rela- 
tivity which allows non-symmetric fields. The theory without a cosmological 
constant was first proposed by Einstein and Straus [iHSllllll] . Schrodinger 
later showed that it could be derived from a very simple Lagrangian density 
if a cosmological constant was included [6^, 7, 8J. Einstein and Schrodinger sus- 
pected that the theory might include electrodynamics, but no Lorentz force 
was found [9lfT0] when using the Einstein-Infeld-Hoffmann (EIH) method pTl 
[T2] . Here we show that a simple modification of the Einstein-Schrodinger the- 
ory closely approximates Einstein-Maxwell theory, and the Lorentz force does 
result from the EIH method, and in fact the ordinary Lorentz force equation 
results when sources are included. The modification is the addition of a sec- 
ond cosmological term A^g^^i,, where f/^t, is the symmetric metric. We assume 
this term nearly cancels Schrodinger's "bare" cosmological term ylfcA^^i,, where 
Nfj^i, is the nonsymmetric fundamental tensor. The total cosmological constant 
A = Ai, + Az can then match cosmological measurements of the accelerating 
universe. Our theory is related to one in [T3], but it is roughly the electromag- 
netic dual of that theory, and it allows coupling to additional fields (sources), 
and it allows Ay^O. 

The origin of our A^ is unknown. One possibility is that A^ could arise from 
vacuum fluctuations, an idea that has been discussed by many authors [M l flS l 
[TBlfTT] . Zero-point fluctuations are essential to both QED and the Standard- 
Model, and are the cause of the Casimir force[T5j and other effects. Another 
possibility is that Az arises dynamically, related to the minimum of a potential 
of some additional field in the theory. Speculation about the origin of this 
second cosmological constant is beyond the scope of this paper. Our main goal 
here is to demonstrate that the theory closely approximates Einstein-Maxwell 
theory. 

Like Einstein-Maxwell theory, our theory can be coupled to additional fields 
using a symmetric metric and vector potential A^ , and it is invariant under 
a U{1) gauge transformation. The theory does not enlarge the invariance group. 
When coupled to the Standard Model, the combined Lagrangian is invariant 
under the usual U{1) (S> SU{2) (§> SU{3) gauge group. The usual U{1) gauge term 
F^'^F^^ is incorporated together with the geometry, and is not explicitly in the 
Lagrangian. Whether this is a step backwards from Einstein-Maxwell theory 
coupled to the Standard Model, or whether the SU{2) and SU{3) gauge terms 
could also be incorporated using non-Abelian fields as in p!?[|T^ . or by using 
higher space-time dimensions, is speculation beyond the scope of this paper. 

This paper is organized as follows. In fj2]we discuss the Lagrangian den- 
sity. In ^j3]-fj5]we derive the field equations and quantify how closely they ap- 
proximate the field equations of Einstein-Maxwell theory. In ^ we derive the 
ordinary Lorentz force equation by taking the divergence of the Einstein equa- 
tions when sources are included. In fj7] we derive the Lorentz force using the 
EIH method, which requires no sources in the Lagrangian. In Sj8] we give an 
exact charged solution and show that it closely approximates the Reissner- 
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Nordstrom solution. In S}9]we give an exact electromagnetic plane- wave solu- 
tion which is identical to its counterpart in Einstein-Maxwell theory. 



2 The Lagrangian density 

Einstein-Maxwell theory can be derived from a Palatini Lagrangian density, 

^^[a,p].9"V^^[,..H • • •)■ (1) 



IGtt 
47r' 



Here Ab is a bare cosmological constant. The Cm term couples the metric g^j/ 
and electromagnetic potential to additional fields, such as a hydrodynamic 
velocity vector m"^, spin-1/2 wavefunction ip, or perhaps the other fields of 
the Standard Model. The original Einstein-Schrodinger theory allows a non- 
symmetric N^i, and in place of the symmetric gfj^i, and , and excludes 
the ^^—gAya.p\9°'^g'^'^A^tL,v\ term. Our "yl-renormalized" Einstein-Schrodinger 
theory introduces an additional cosmological term ^—gAz , 



iDTT L 

where A^f^—A^ so that the total A matches astronomical measurements . 

A^Ab + AzP^ IQ-^^cm-^, (3) 
and the physical metric and electromagnetic potential are defined to be 



iOTT 



A,. 



r[:,]/[(n-l)v^2^]. 



(4) 



Eq. dH) defines g'^" unambiguously because ~ [—det{ 

Here and throughout this paper we use geometrized units with c — G — 1, 
the symbols ( ) and [ ] around indices indicate symmetrization and antisym- 
metrization, g = det{gp_,y), N = det{N^^), and N^'^^ is the inverse of N^^ such 
that N^'^^Nyp^ = 5'^ . The dimension is assumed to be n=4, but "n" is retained in 
the equations to show how easily the theory can be generalized. The Cm term 
is not to include a \/—gA[a,p]g°''^g^'^A[fj^ ,j] part but may contain the rest of the 
Standard Model. In ([2]), TZu^iT) is a form of Hermitianized Ricci tensor[T], 



^I^H^Za) ^ r^a^afi^ -^riy] -^q^] /("-^ " ("'^ 

0, 



= Oandr„"[,_^j: 



This tensor reduces to the ordinary Ricci tensor when 
as occurs in ordinary general relativity. 

It is helpful to decompose F"^ into a new connection F"^ , and A^r from ([4]) 



where F^^^FZ^ 



(6) 
(7) 
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By contracting ([7|) on the right and left we see that -T"^ has the symmetry 
SO it has only rt' — n independent components. Substituting ([6]) into ([5]) gives 



7^,^(^) = 7^,^(f) + 2^[,,^]^2^. (9) 

Using ([9]), the Lagrangian density ([2]) can be written in terms of F'^^ and A^, 



C{rJ'r,Npr) =-^^ [iV^'"'(7t.^ + 2^[,,^]v/^) + {n-2)Ab 

-T^ {n-2)A, + 7/., g^,, A,...). (10) 

Here iZ^f_i=TZ^^{r), and from (|8l5p we have 

From (|6|8p . /""^ and fully parameterize F"^ and can be treated as in- 
dependent variables. It is simpler to calculate the field equations by setting 
6£/SF^^ = and dC/SA,^ — instead of setting 5C/SF^^ = 0, so we will follow 
this method. 

To do quantitative comparisons of this theory to Einstein-Maxwell theory 
we will need to use some value for Az. One possibility is that Az results from 
zero-point fluctuations! 1 4 l [T5 |l 1 6U 1 7] . in which case using ([3]) we get 

Ab « ~Az ~ Czojjll - lO^'^cm-^ (12) 

ujc = (cutoff frequency) ^ 1/Zp, (13) 

(J — _L ( fcrmion _ boson \ ^ 60 
^ ^ 2tt Vspin states spin states/ 2tt 

where Ip = (Planck length) ~ 1.6 x lO^'^^cm. We will also consider the limit 
ujc—^oo, \Az\— I- oo, ylb^oo as in QED, and we will prove that 

lim / yl-renormalized \ _ /Einstein-MaxwellX /-, rN 

jylzl^cx) V Einstein-Schrodinger theory / ^ V theory /' ^ ' 

The Hermitianized Ricci tensor ^ has the following invariance properties 

n^^,{r^) = TZf^^iF) (T = transpose), (16) 

n,^if^,+ Sfp^,r]) = 7^,^(^p" ) for an arbitrary ^{x^). (17) 
From (|16ll7p . the Lagrangians (|2ll0p are invariant under charge conjugation, 
Q^-Q, A„^-A,, r^^^r^,, Cp^/;"., N.^^N,,, iV^''^^iV^^^(18) 
and also under an electromagnetic gauge transformation 

2n 



i^^^e't A^^A^-^^,^, r^,^r^,, r;^-^r;,+^<5p^0,,]v^2^, (19) 
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assuming that Cm is invariant. With yl;, > 0, yl^ < as in (fT2|) then J""^, r""^, 

A^jy^ and N^'^^ are all Hermitian, TJ-jy^ and TZ^^{r) are Hermitian from (I16p . 
and (?jy^, ^o- and £ are real from (|4l2ll0p . 

In this theory the metric ^ is used for measuring space-time intervals, for 
calculating geodesies, and for raising and lowering of indices. The covariant 
derivative ";" is always done using the Christoffel connection formed from g^^, 

Cm = \ 5""(5m-.- + 5^-.^^ - (20) 

We will see that taking the divergence of the Einstein equations using (|20l4p 
gives the ordinary Lorentz force equation. The electromagnetic field is defined 
in terms of the potential (j!]) 

F^v = Ay^^ — Afj^^u. (21) 
However, we will also define another field /'"^ 

y/^f"" - V^N^^'^^'^Al^^/v^i. (22) 
Then from g]), g''" and f'"'y/2iA^^^^'^ are parts of a total field, 

{V^/^)N^'''' = g^'' + /'"'%/2a;^/l (23) 

We will see that the field equations require fp,i,~F^^ to a very high precision. 
The definitions ^ of and A^ in terms of the "fundamental" fields Npr , 
may seem unnatural from an empirical viewpoint. On the other hand, our 
Lagrangian density ([2]) seems simpler than ([T]) of Einstein-Maxwell theory, it 
contains fewer fields, and these fields have no symmetry restrictions. However, 
these are all very subjective considerations. It is much more important that 
our theory closely matches Einstein-Maxwell theory, and hence measurement. 

Note that there are many nonsymmetric generalizations of the Ricci tensor 
besides the Hermitianized Ricci tensor TZyf_i{F) from ^ and the ordinary 
Ricci tensor R^p{r). For example, we could form any weighted average of 
Riyf_i[F), Rf^^{r), R^^{r^) and Rfj.p{r^), and then add any linear combination 
of the tensors f " ,, and T,°' .f," All of these 

generalized Ricci tensors would be linear in F^^^, quadratic in T"^, and would 
reduce to the ordinary Ricci tensor when = and F2\^^ ^ occurs in 

ordinary general relativity. Even if we limit the tensor to only four terms, there 
are still eight possibilities. We assert that invariance properties like (|16ll7p are 
the most sensible way to choose among the different alternatives, not criteria 
such as the number of terms in the expression. 

Finally, let us discuss some notation issues. We use the symbol /^'^ for the 
Christoffel connection ([20]) whereas Einstein and Schrodinger used it for our 
-T"^ and -Z^" respectively. We use the symbol g^v for the symmetric metric 
(HI) whereas Einstein and Schrodinger used it for our A^^^, the nonsymmetric 
fundamental tensor. Also, to represent the inverse of N^fj. we use 7V~'°'" instead 
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of the more conventional A''"'^, because this latter notation would be ambigu- 
ous when using g^^'^ to raise indices. While our notation differs from previous 
literature on the Einstein-Schrodinger theory, this change is required by our 
explicit metric definition, and it is necessary to be consistent with the much 
larger body of literature on Einstein-Maxwell theory. 



3 The Einstein equations 



To set SC/S{\/—NN^'^'') = we need some initial results. Using Q and the 
identities det{sM) = s"det{M), det{M-^)^ l/det{M) gives 

y/^ = {-det{V^N^^"')y/^''-^\ (24) 
= (-det(V^5^''))i/("-^) = {-det{^/^N^'•^"'^))^/'•'''^\ (25) 

Using ()24l25l4|) and the identity d{det{M-)) / dAP'" ^ M-^det{M-) gives 



N, 



Setting SC/6{y/^N~^'"')= using (|10l26p gives the field equations. 



where Si^fj, and the energy-momentum tensor Ti,^ are defined by 



Si/ 



1 



T = ^ — — n <?" 



C — T n T" 



(n-2)- 



(26) 

(27) 

(28) 
(29) 



The second equality in (|28p results because £,„ in ([2]) contains only the metric 
y^gt^'' = y/^N~^^'"'^ from (gl), and not V^N'^^t^'^l Taking the symmetric 
and antisymmetric parts of ([27]) and using ([2T|) gives 



J^^f^.T:) , (30) 
(31) 



(32) 



Also from the curl of (|3ip we get 



To put (j30p into a form which looks more like the ordinary Einstein equa- 
tions, we need some preliminary results. The definitions (|4l22p of g^^ and fi,^ 
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can be inverted exactly to give N^^ in terms of and f^^. An expansion in 
powers of A^^ will better serve our purposes, and is derived in Appendix [Aj 

iV(.M)=.g.M - 2(^/."/.M - ^^(^)9.,ru}jA-^,^ + {f)A-^ . . . (33) 

N^u^.]=U^.V2^Al^'^ + (/3)yi;=^/' .... (34) 

Here the notation {p) and (/*) is for terms like fi^af^af^tj. and fuaf^crf^pffj.- 
Let us consider the size of these higher order terms relative to the leading order 
term for worst-case fields accessible to measurement. In geometrized units an 
elementary charge has 

Q^ = eJ^■^ = Valp = 1.38 x IQ-^^cm (35) 

V \ nc c'^ 

where a — e'^ /he is the fine structure constant and Ip = ^ Gfij (? is the Planck 
length. If we assume that charged particles retain /^o ~ Q/?"^ down to the 
smallest radii probed by high energy particle physics experiments (10~^^cm) 
we have from (|35ll2p . 

l/'oPMfc - (Qe/(10-i')2)VA ~ 10-s^ (36) 

Here |/^o| is assumed to be in some standard spherical or cartesian coordi- 
nate system. If an equation has a tensor term which can be neglected in one 
coordinate system, it can be neglected in any coordinate system, so it is only 
necessary to prove it in one coordinate system. The fields at 10~^^cm from an 
elementary charge would be larger than near any macroscopic charged object, 
and would also be larger than the strongest plane-wave fields. Therefore the 
higher order terms in (|33ll34p must be < 10^®^ of the leading order terms, so 
they will be completely negligible for most purposes. 

In fj5]we will calculate the connection equations resulting from SC/SF^^ — Q. 
Solving these equations gives (|53I54I57I59|) . which can be abbreviated as 

^r.,) = Cm + 0{A-,'). Tf^^j - 0{A-"^), (37) 
= 0{A-,^), 7^[,^] = 0{A-^'^), (38) 

where F^^ is the Christoffel connection ([20)l . TZi,i_, =TZijf_,{r), Rvfi^ Rvti[r) and 

= - ^ gutJ.Tip, Gufj. = Riyf_i - | gu^iR- (39) 

In dnSD the notation 0(yl^^) and 0(yl^^^^) indicates terms like f ^-af f,;aA^^ 
and /[^M,a];"^b^^^- 

From the antisymmetric part of the field equations ([3T|l and (|34|38p we get 



U = F,^ + 0{Al'). 



(40) 
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So fi,^ and F^^ only differ by terms with Ab in the denominator, and the two 
become identical in the limit as Ab—^oo. Combining (|30p with its contraction, 
and substituting p9|33|3p gives the Einstein equations 

= SttT,^ - At (^iV(,^) - ^g.^iV;^ +A,(^~l^g,^, (41) 
= 8nT,^ + 2 (^Ul.^-\9.^ru)j + ^ (f - l) 9^, + • • • -(42) 

From (|28|29|) we see that T^^ will be the same as in ordinary general relativ- 
ity, for example when we include classical hydrodynamics or spin- 1/2 fields 
as in [2T1l22j. Therefore from (|40|38p . equation ((42|) differs from the ordinary 
Einstein equations only by terms with Af, in the denominator, and it becomes 
identical to the ordinary Einstein equations in the limit as ylb — *■ oo (with an 
observationally valid total A). In ^we will examine how close the approxi- 
mation is for Ab from (fT2|) . 



4 Maxwell's equations 



Setting 5C/6Ar = and using (|10|22p gives 


where 

i 



-1 



-9 



dCm ( dCr. 



dAr V^^r, 



From (|43l2ip we get Maxwell's equations. 



(43) 



(44) 



(45) 
(46) 



where /^.^ = F^^, + 0(yl^^) from (gO]). From ((2]44)) we see that f will be the 
same as in ordinary general relativity, for example when we include classical 
hydrodynamics or spin- 1/2 fields as in pTll^ . From we see that equations 
(|45l46p differ from the ordinary Maxwell equations only by terms with A^ 
in the denominator, and these equations become identical to the ordinary 
Maxwell equations in the limit as A^^ oo. In fj5] we will examine how close 
the approximation is for A}, from (|12p . 

Because Lm in © couples to additional fields only through g^i, and A^, 
any equations associated with additional fields will be the same as in ordinary 
general relativity. For example in the spin- 1/2 case, setting 6£/64> = will 
give the ordinary Dirac equation in curved space as in pTU^ . It would be 
interesting to investigate what results if one includes /^^ , N^^^ 



or in 
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although there does not appear to be any empirical reason for doing so. A 
continuity equation follows from (|45|) regardless of the type of source, 

/;p-^r^[r;.]=0. (47) 

Note that the covariant derivative in (|45l47p is done using the Christoffel 
connection (j20p formed from the symmetric metric 



5 The connection equations 



Setting S£/6r^^ — with a Lagrange mulitiplier term to enforce the 

symmetry ([5]), and using (|10l43p gives 



,1/2 



, ^, .... -^^'"S^- (48) 

These are the connection equations, analogous to g^"^ -.p = in the symmet- 
ric case. Note that we can also derive Ampere's law (|43p by antisymmetriz- 
ing and contracting these equations. From the definition of matrix inverse, 
N^p-'={\/N)dN/dNrp and N^P^Nr^, = 5P we get the identity 



-TV' 



dNrp 



-N 



N^P\aNrp. (49) 



Contracting ([48]) with N^p using (|8l49p . and dividing this by (n — 2) gives. 



(V^),^-C^^: 



{n-l)[n-2)Al'^ 



(50) 



Multiplying (gS]) by -N^pNrp and using ^ gives 



Nvp.,13 - r^^pNa^ - FppN^o 



N[aP]Ni,fj\ . 



(n-2) 



r- (51) 



Equation ([5T|) together with (|30l32l8p are often used to define the Einstein- 
Schrodinger theory, particularly when T,yp — 0, j" = 0. 

Equations or ([5T|) can be solved exactly similar to the way 

9pT-p — can be solved to get the Christoffel connection. An expansion in 
powers of A^^ will better serve our purposes, and such an expansion is derived 
in Appendix E of [H], and is also stated without derivation in [25] . 

(52) 



^[up]- 



p)+j^^{{rup),''9up-2{rup).^j';,^) 



47r 



[n-2y 

Tiifi^p; +/ p-.f^f v:p) + p] 



v + (53) 



^/2^A-"^+{f')A~"\.{^A) 
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In JSg), is the Christoffel connection The notation (f) and (J^') 

refers to terms hke rrFafM and frF pfy.-t.)- As in pl34D . we 
see from (pB)) that the higher order terms in (|53ll54p must be < 10^^^ of the 
leading order terms, so they wiU be completely negligible for most purposes. 
Extracting TJ^ of from the Hermitianized Ricci tensor pip gives, 

'<-(i^m)1^ -n-iy^i^ J + -'(l/^);Q~-'Q(iy;^)~-'(iyQ)-'(o-^)~-'[iya]-'[CT/j]+-'(l/M)-'(TQ' lOOj 

l<-[vti]y^ )— J-liyfi]-a^-'-{iya)-'-[afi]^-'-liya]-'-(afj.)+-'-[^fj.]-'-cra- 

Substituting (|52j|54l45]) into ([55| . and using ([351) gives 
1 s 

fTfS r a ( fpar \a ^ f f[<^P o<\ 

iJupJ JP ■,T-,a ^yup.\J Jap), ;0! ^yi^pJ[af3;a]J ; 

+ 8^J - (f^J-^M + ^^^9^,fjp + (/')) • • ■ ■ (57) 

From (|42p we can define an "effective" energy momentum tensor T^^^ which 
applies when G^^ is used in the Einstein equations and Cm — 0, 

Snf,^ = 2 (^/.7.M- J 9u^np^ - (G.^ - G,^). (58) 
Substituting (|54l45p into dM]) gives 

''^M = ^2'^['''''"1^"^^-^"^'['''"1~^-^"'''['''"1 ■ ■ ■ (59) 

As we have noted in Sj3]and ^ the ylb in the denominator of (|57l59p causes 
our Einstein and Maxwell equations (|42l45l46p to become the ordinary Einstein 
and Maxwell equations in the limit as o^c^oo, \Az\—^oo, A{,—^c>o, and it also 
causes the relation fyp,~Fy^ from to become exact in this limit. Let us 
examine how close these approximations are when Af, ~ 10^^cto~^ as in (|12p . 

We will start with the Einstein equations ((42|) . Let us consider worst-case 
values of the 0{A^^) terms in ([57)) and compare these to the ordinary electro- 
magnetic term in ([58)1 . If we assume that charged particles retain f^o^Q/r^ 
down to the smallest radii probed by high energy particle physics experiments 
(10~^^cm) we have, 

Ifoa/fWVA, ^ 4Mb (10-1^)2 ^ 10-32, (60) 
\foaa/f\\/Ab ^ 6/ A, {lO-^'f ^ 10^^^ (61) 

So for electric monopole fields, terms like J" v.af" p.;<tA^'^ and f^^ fr(i^: p);aA^^ 
in ([57)1 must be < lO^^^ of the ordinary electromagnetic term in ([55)) . And 
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regarding j"^ as a substitute for (l/47r)/"'^.^ from (|45|) . the same is true for 
the j"^ terms. For an electromagnetic plane-wave in a flat background space 

^ Ae^sm{kaX°') , e"e„ = -1 , A:"fc„ = k^e^ = 0, (62) 
/.M = 2A[^^,] = 2Ae[^A:,]Cos(fc„a;"), f = 0. (63) 

Here A is the magnitude, k" is the wavevector, and is the polarization. 
Substituting (|62|63p into ([57|) . all of the terms vanish for a flat background 
space. Also, for the highest energy gamma rays known in nature (10^'^eV, 
lO^^Hz) we have from (HH), 

\f\;i/.f\\VA, ^ {Elficf/A, ^ 10-l^ (64) 
\f\;i;i/f\\IA, ^ {E/hcf/A, ^ 10-16. (65) 



So for electromagnetic plane- wave fields, even if some of the terms in ((5^ 
were non-zero because of spatial curvatures, they must still be < 10"^^ of the 
ordinary electromagnetic term in (I58|) . Therefore even for the most extreme 
worst-case fields accessible to measurement, the extra terms in the Einstein 
equations (|42|) must all be <\Q~^^ of the ordinary electromagnetic term. 

Now let us look at the approximation f^^ w F^^ from , and Maxwell's 
equations (|45I46|) . From the covariant derivative commutation rule, the cyclic 
identity 2_Ry[^ct]p — Ri^fiar, the definition of the Weyl tensor C^^aT, and the 
Einstein equations Ry^^ — —Ag^f^ + {p) . . . from ((42|) we get 

2/"zy;[p;Q] = -^f^^CaTiyfl, + /t^M + (/^) (66) 

Substituting ([M)) into the field equations ([?T|) and using (|59l66p we get 



(n-1) 

where Etu^iol — (Levi— Civita tensor), Carvn — (Weyl tensor), and 
1 2 

— '^f[ufi,a]^T ^ I f[iyfi.a] ~ 3 vfJ.a- (68) 

The 9[T.a]^i^ti'^"A^^ term in ([57]) is divergenceless so that it has no effect on 
Ampere's law gt]). The U^^A/Ab term is ~ 10-^22 fj.^^^ (j3TT2)) . The 

(/■^)yl^i term is < 10"^^ of /j,^ from ([36]) . The largest observable values of 
the Weyl tensor might be expected to occur near the Schwarzschild radius, 
Vs — 2Gm /c^, of black holes , where it takes on values around Vs/r^ ■ The largest 
value of rg/r^ would occur near the lightest black holes, which would be of 
about one solar mass, where from (|12p . 



Ab ^r2 Ab \2Gm ' " ^ ' 
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And regarding j"^ as a substitute for (l/47r)/'^'^.^ from ((45|) . the term 
is < 10"'^^ of f^fj^ from ((6T|) . Therefore, the last four terms in ([67| must all 
be < 10^'^^ of /y^. Consequently, even for the most extreme worst-case fields 
accessible to measurement, the extra terms in Maxwell's equations (|45l46p 
must be < 10~^^ of the ordinary terms. 

The divergenceless term 9\^^^a-\eufJ"A^^ of (p7|) should also be expected to 
be < 10~^^ of /^^ from (|60l61l68p . However, we need to consider the possibility 
where 9t changes extremely rapidly. Taking the curl of ((67)) . the and 
terms drop out, and we get something similar to the Proca equation|26p27j . 

Op = (- e,, . - + ie/'^^(/"^C„.[.,)^.] + ^^^^ dp + (/'')) ^ • • ■ ■ (70) 

Here the constraint = results from the definition ([68]) and we are using 
a (1,-1,— 1,-1) metric signature. Eq. ifTH)) suggests that Op Proca-wave so- 
lutions might exist in this theory. Assuming that the magnitude of Carup, is 
roughly proportional to Op for such waves, and assuming that /^^ goes accord- 
ing to ([ST]) with Ffj^^ — O, the extra terms in ([70]) could perhaps be neglected in 
the weak field approximation. Using ([70)1 and Ai,Ki —A^ = CzUJ'^lp from (|12p . 
such Proca-wave solutions would have an extremely high minimum frequency 

wproca = \/2^ «\/2C^w2^P ~ IQ-^Vad/s, (71) 

where the cutoff frequency uJc and Cz come from (|13ll4p . 

There are several points to make about (|70l7ip . 1) A particle associated 
with a Op field would have mass hujproca, which is much greater than could be 
produced by particle accelerators, and so it would presumably not conflict with 
high energy physics experiments. 2) We have recently shown that sin[kr^Lut\ 
Proca-wave solutions do not exist in the theory, using an asyptotically flat 
Newman- Penrose l/r expansion similar to [28ll29j. However, it is still possible 
that wave-packet solutions could exist. 3) Substituting the k — flat space 
Proca-wave solution Op = (0, 1,0, O)sin[ujprocat] and Fp^ — into (|67l58l57p . 
and assuming a flat background space gives Too = — 2/ylf, < 0. This suggests 
that Proca-wave solutions might have negative energy, but because sin[ kr—ojt] 
solutions do not exist, and because of the other approximations used, this 
calculation is extremely uncertain. 4) With a cutoff frequency LOc^l/lp from 
(fT3| we have ojproca > from (|71|13|14p . so Proca-waves would presumably 
be cut off. More precisely, ([7T|) says that Proca-waves would be cut off if 
ojc > 1/(^pV2C'z ). Whether uJc is caused by a discreteness, uncertainty or 
foaminess of spacetime near the Planck length[30,31,32,33,34j, or by some 
other effect, the same uJc which cuts off Az in ([T2|) should also cut off very 
high frequency electromagnetic and gravitational waves, and Proca-waves. 5) 
If wave-packet Proca-wave solutions do exist, and they have negative energy, 
it is possible that Op could function as a kind of built-in Pauli-Villars field. 
Pauli-Villars regularization in quantum electrodynamics requires a negative 
energy Proca field with a mass hujproca that goes to infinity as — > oo, as 
we have from (|7ip . 6) As mentioned initially, it might be more correct to 
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take the limit of this theory as Wc — > oo, jylzj — > oo, ylb — s- oo, as in quantum 
electrodynamics. In this limit (|70|7ip require that 0p— >0 or uoproca^oo, and 
the theory becomes exactly Einstein-Maxwell theory as in p5|l . 7) Finally, we 
should emphasize that Proca-wave solutions are only a possibility suggested 
by equation ((70)) . Their existence and their possible interpretation are just 
speculation at this point. We are continuing to pursue these questions. 



6 The Lorentz force equation 

A generalized contracted Bianchi identity for this theory can be derived using 
only the connection equations (|48| and the symmetry (jS]) of J""^, 

{^PnN^'""R„x + ^FNN^'""iZxa)M - ^^NN^^^Hau^x = 0. (72) 
This identity can also be written in terms of g'*'^, fP'^ and G^^ from (I4l22l39p . 

GZ;, = {lr'%p,.]+n;^%.])V2iA-'^\ (73) 

The identity was originally derived [3l[7] assuming j'^ = in (|48p . The derivation 
for j'^T^O was first done[25j by applying an infinitesimal coordinate transfor- 
mation to an invariant integral, and it is also done in Appendix B of ^21j using a 
much different direct computation method. Clearly ()72I73|) are generalizations 
of the ordinary contracted Bianchi identity 2{y/~g W — \/—g g^'^Rav,\ =0 
or GJ^.g-^O, which is also valid in this theory. 

Another useful identity[T3] can be derived using only the definitions (|4l22p 

(N^\^-h>>tNi),^ ^ (f r^^^kP.H + ^',aN^,.^^^A-,^'\ (74) 

The ordinary Lorentz force equation results from taking the divergence of 
the Einstein equations jH]) using (|73|45|31|74|2ip 



rpa ^ 

'^''^ ~87r 



yl,(7V(%)-|<5^;A^;),J (75) 



Note that the covariant derivatives in (|73l74l75p are all done using the Christof- 
fel connection (|20p formed from the symmetric metric ([4]). 



7 The Einstein-Infeld-Hoffmann Equations of motion 

For Einstein-Maxwell theory, the EIH method allows the equations of mo- 
tion to be derived directly from the electro-vac field equations. For neutral 
particles the method has been verified to Post-Newtonian order [TT], and in 
fact it was the method first used to derive the Post-Newtonian equations 
of motion p5]. For charged particles the method has been verified to Post- 
Coulombian order[12,36,37j, meaning that it gives the same result as the Dar- 
win Lagrangian|27j. In we derived the exact Lorentz force equation for 
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this theory by including source terms in the Lagrangian. Here we derive the 
Lorentz force using the EIH method because it requires no source terms, and 
also to show definitely that the well known negative result of [QllTO] for the 
unmodified Einstein-Schrodinger theory does not apply to the present theory. 
We will only cover the bare essentials of the EIH method which are necessary 
to derive the Lorentz force. We will also only calculate the equations of motion 
to Newtonian/ Coulombian order, because this is the order where the Lorentz 
force first appears. 

The EIH method assumes the "slow motion approximation" , meaning that 
v/c^l. The fields are expanded in the form[TT | fT2 ( l36 p 37]. 



9tJ.u 




(76) 


700 


— 2700^^^ + 4700^^^ . . . 


(77) 


7ofc 


— SlOk^^ + 570fcA^ . . . 


(78) 


lik 


= 47ifeA . . . 


(79) 




= + 4A0A4 . . . 


(80) 


Ak 


= s^fcA^ + s^feA^ . . . 


(81) 




= 2/ofcA^ + 4/ofeA^ ■ • ■ 


(82) 


fik 


— sfik^^ + 5/ifcA^ . . . 


(83) 



where A ~ w/c is the expansion parameter, the order of each term is indicated 
with a left subscript '9] , ?7^j/ = diag(l, —1, —1, —1), and Latin indices run from 
1-3. The field 7^^ (often called h^^ in other contexts) is used instead of 5^1/ only 
because it simplifies the calculations. Because \ ^ v/c, when the expansions 
are substituted into the Einstein and Maxwell equations, a time derivative 
counts the same as one higher order in A. The general procedure is to substitute 
the expansions, and solve the resulting field equations order by order in A, 
continuing to higher orders until a desired level of accuracy is achieved. At 
each order in A, one of the ;7^j, terms and one of the if^n terms will be 
unknowns, and the equations will involve known results from previous orders 
because of the nonlinearity of the Einstein equations. 

The expansions (|77II83I) use only alternate powers of A essentially because 
the Einstein and Maxwell equations are second order differential equations |35j. 
although for higher powers of A, all terms must be included to predict radiation|12t 
136.37 . Because A~w/c, the expansions have the magnetic components Ak and 
fik due to motion at one order higher in A than the electric components Ao 
and /oi- As in [T2l[36ll37] . fak and fik have even and odd powers of A respec- 
tively. This is the opposite of 9,10 because we are assuming a direct defi- 
nition of the electromagnetic field p2|34|67|21|) instead of the dual definition 
jap ^ £"P'^M^rj^^j/2 assumed in [9l[T0|. 

The field equations are assumed to be of the standard form 

G^. = 87rT,,, where G^. = i?,,. - ig^.g^'^i^a/s, (84) 
or Rf^^ = %TTSf^ty where S^,y ^ T^^, ~ ^g^jyg"'^ T^p- (85) 
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However, with the EIH method we solve a sort of quasi-Einstein equations, 

= G^. - 87rf^„ (86) 

where 

Gp,i, — R^y — ■^rjf^^ri"'^ Raf3, T^y — S^y — ^rif^yT]"'^ Saf3- (87) 

Here the use of rj^i, instead of g^i, is not an approximation because (j85p impUes 
(f86l) whether G^i/ and T^i, are defined with -q^i, or g^i/. Note that the references 
use many different notations in (|86|) : instead of G^i^ others use 77^^/2 + yl^i^, 
^nuf^+A^y or [LS:^i^] and instead of SnTf^y others use —25*^^, —A'^y, —A^^^ 
or [RS : ^iv] . 

The equations of motion result as a condition that the field equations (|86[) 
have acceptable solutions. In the language of the EIH method, acceptable solu- 
tions are those that contain only "pole" terms and no "dipole" terms, and this 
can be viewed as a requirement that the solutions should resemble Reissner- 
Nordstrom solutions asymptotically. To express the condition of solvability 
we must consider the integral of the field equations ([5S)) over 2D surfaces S 
surrounding each singularity. 



1 



s 



iCn = — / { iG^k - S-K iTi_,k)nkdS. 
Ztt J 

Here is the surface normal and I is the order in A. Assuming that the 
divergence of the Einstein equations (|84p vanishes, and that (|86|) has been 
solved to all previous orders, it can be shownjTI] that in the current order 

{iG^k~8niff,k)\k = 0- (89) 

Here and throughout this section "|" represents ordinary derivative TT. From 
Green's theorem, ((5^ implies that in ([55)1 will be independent of surface 
size and shapeHTl. The condition for the existence of an acceptable solution 
for 47,jt is simply 

iC, = 0, (90) 

and these are also our three 0{X'^) equations of motion|ll|. The Co component 
of (|88p causes no constraint on the motion [TT] so we only need to calculate 
Gik and T^fe. 

At this point let us introduce a Lemma from llj which is derived from 
Stokes's theorem. This Lemma states that 

^{■■■)ki\inkdS = Q if T(...)ki^ ~^(-)ik, (91) 

where !F(...)m is any antisymmetric function of the coordinates, Uk is the surface 
normal, and S is any closed 2D surface which may surround a singularity. The 
equation 46*,; = is a condition for the existence of a solution for 47^^ because 
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47ife is found by solving the O(A^) field equations (|86|l . and is the integral 
([88]) of these equations. However, because of the Lemma (|9T|) . it happens 
that the 47ife terms in 4Gik integrate to zero in ([88| . so that 4(7^ is actually 
independent of 47^^. In fact it is a general rule that Ci for one order can be 
calculated using only results from previous orders [Tl]. and this is a crucial 
aspect of the EIH method. Therefore, the calculation of the O(A^) equations 
of motion (|90p does not involve the calculation of 47^^, and we will see below 
that it also does not involve the calculation of 3/^^ or 4/ofc. 

The iGik contribution to (l88|) is derived in llj. For two particles with 
masses nii, TO2 and positions <^2, the 0{\*) term from the integral over the 
first particle is 

'tC- = ^j iG,knkdS = -4 jmili - mims^ | , (92) 

where 



^ = VW^WK^)- (93) 

If there is no other contribution to (j88p . then (j90p requires that'fa=0in ((92|) . 
and the particle acceleration will be proportional to a V(mim2/r) Newtonian 
gravitational force. These are the EIH equations of motion for vacuum general 
relativity to ©(A^), or Newtonian order. 

Because our effective energy momentum tensor (|58p is quadratic in /^^, and 
the expansions ()77ll83p begin with A^ terms, the 0{\^) — 0{\^) calculations 
leading to ([92]) are unaffected by the addition of the electromagnetic terms to 
the vacuum field equations. However, the 8t: ^Tik contribution to (j88p will add 
to the ^Gik contribution. To calculate this contribution, we will assume that 
our singularities in f^^ are simple moving Coulomb potentials, and that 0^ = 0, 
A^O. Then from (|67l82j|85|) we see that 2^ofc —2foki a-nd from inspection of the 
extra terms in our Maxwell equations (|45l46l67p and Proca equation (ffO)) . we 
see that these equations are both solved to 0{\'^). Because ([551) is quadratic 
in /^i,, we see from (j82ll83l) that only 2/ofc can affect the 0{X^) equations of 
motion. Including only 2/ofci our /^^ is then a sum of two Coulomb potentials 
with charges Qi, Q2 and positions Cij ^2 of the form 

2^^ = (2¥', 0,0,0) , 2/0fc=2 2A[fc|0] =-2V|fc, (94) 

2(p = V''+V' , ^^^Qi/ri , i/''=Q2/r2, (95) 
ra= V(^^-C)(^^-C) ' a-1...2. (96) 

Because ([58| is quadratic in both ffj^^ and 5^,^, and the expansions ([771 - 
start at A^ in both of these quantitites, no gravitational-electromagnetic 
interactions will occur at 0{X'^). This allows us to replace covariant derivatives 
with ordinary derivatives, and with r/^^ in ([55]). This also allows us to 
replace T^^ from (|86l87p with ((SS]). Keeping only 0{X*) terms when ([M]) is 
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substituted, the spacial part of ([58| gives, 

/o(5|m)|a + / s|a/om| +/ a|s/ 0|m 2 '?^™ 

Note that 2<y3 from ([55)1 obeys Gauss's law, 

V\a\a^O. (98) 

Substituting ^ into ((971) and using ^ gives 



(99) 



87r4f;„i = -2 ^<y3|s(y5|„i + -?7sm'^|r<y5|r 

s|m|a '/'|s|a'/'|m|a I VlalsVlalm I r) '^swi l¥'|rV-'|r j |a|a I ^^6 

From (PT|) . the second group of terms in pOOp integrates to zero in so 
it can have no effect on the equations of motion. The first group of terms in 
(|100P is what one gets with Einstein-MaxweU theorv [T2l[36ll37j . so at this stage 
we have effectively proven that the theory predicts a Lorentz force. 

For completeness we will finish the derivation. First, we see from (|100l98p 
that 4Tg„-^^s = 0. This is to be expected because of and it means that 

the Sir 4Tsm contribution to the surface integral (|88|) will be independent of 
surface size and shape. This also means that only 1/distance^ terms such as 
'Hsm/f'^ or XsXm/'r'^ can contribute to (|88p . The integral over a term with any 
other distance-dependence would depend on the surface radius, and therefore 
we know beforehand that it must vanish or cancel with other similar terms 1 1]. 
Now, = from ((95)) . Because ip^^ and 'ip'\i both go as 1/distance^, but 

are in different locations, it is clear from (|100p that contributions can only 
come from cross terms between the two. Including only these terms gives, 

8^4^™ = -2(^|.V'|„. + ^|.V'|„.+ Vsmi^\r^l). (101) 



Some integrals we will need can be found in [llj . With ij] — 1/ y/x^liF we have, 

1 /-V 1 rV ... 1 



, , i^lmnmdS = -1 , — / IplaU.mdS = --5am- (102) 

An J ' An J ' 6 

Using (|101ll02l95p and integrating over the first particle we get, 

^ii^ls^lm+i^ls^lm+ Vsmi^lri^frhmdS (103) 

4QiV^f.(a)(-|-l + |) =-4QiV^2(^^)_ (^04) 



1 

2^ 









/■ 


— 87rrs„i 


nmdS=^j 
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Using (|90|88|104I92I95P we get 

= 4^ = -4 jmill - mim2^ (^^^ | - ^Qi^J^^i) (105) 

- "4 jmiCl - [I) + (J) } , (106) 

where 

These are the EIH equations of motion for this theory to 0{X'^), or Newtonian/ 
Coulombian order. These equations clearly exhibit the Lorentz force, and in 
fact they match the C'(A'*) equations of motion of Einstein-Maxwell theory. 



8 An exact electric monopole solution 



Here we give an exact charged solution for this theory which closely approx- 
imates the Reissner-Nordstrom solution [381139] of Einstein-Maxwell theory. A 
MAPLE program[40j which checks the solution and the derivation [41] are 
available. The solution is 



cadt^ dr^ 



1 

ca 



cr'de' 

2 , 



-2-2 

cr stn 



/—N = r sin 9, 



-A' 

— ^2 



1 



AM 4A 



Abr^ SAb 
2M Ar^ Q2y 
3~ ^ 



2 c- 1 - 



AbT^ 



A 



1 



A 
~Ab 



(108) 
(109) 
,(110) 

(111) 



r 6 r 

where (') means d/dr, and c and V are very close to one for ordinary radii 



{2€, 



Ay,r^ 



V 



rAb 



r cdr 



AbT^ 

3 ~ 



[i!]24»(2i-l) V^r4 
n2 



(112) 



1 



{2^) 



10^r4 



and the nonzero connections are 
4a2Q2 



^ 00 



^ 1 ' 



2 Abr^ 

T^2 rn3 7^3 1 

^ 21 — ^ 13 — ^ 31 — „ ' 



pO _ pO 
10 ^ 01 



a 

2^ 



_(2Q 

i!(i + l)!4'(4i + l) \Abr^ 

J ^ 11 



2\ « 



^13) 



—a 



(114) 



fi, = -ar , = -ar sin^fl , fl, = /if, = cot ( 



/"Ij = — sin6'cos6'. 



rn2 



p2 



f3 
30 



aV2iQ -,1 __2aV2iQ 

^ 10 ~ ^ 01 ^ 



br' 
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The solution matches the Reissner-Nordstrom solution except for terms 
which are negligible for ordinary radii. To see this, first recall that A/Ai, ~ 
-[^Q-122 fj-Qjjj (j3|i2p . so the A terms are all extremely tiny. Ignoring the A 
terms and keeping only the 0{A'^^) terms in (|110lllllll2lll3p gives 



Q 



AM 4Q2 



Aq — — 

r 



a = 1 



AhV^ Ahr^ 
M AQ^ 



Abr^ Sylfcr^ 



2M Q2 



IQAbr 



0{A-, 
0{A-, 

+ 0{A-% 



c = 1 - 



+ 0{A-^). 



(115) 
(116) 
(117) 
(118) 



For the smallest radii probed by high-energy particle physics we get from p6p , 
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-66 



(119) 



The worst-case value of M/A\jr^ might be near the Schwarzschild radius 
of black holes where r = rs= 2M and M/Abr^ = l/2Af,r'^. This value wiU be 
largest for the lightest black holes, and the lightest black hole that we can 
expect to observe would be of about one solar mass, where we have 
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M 



Abr^ 2Abrl 2Ab \2GmQ 



'10 



-77 



(120) 



From pi9ll20l3ll2p we see that our electric monopole solution (jlOSIIllip 
has a fractional difference from the Reissner-Nordstrom solution [38l[39] of at 
most 10^®^ for worst-case radii accessible to measurement. Clearly our solu- 
tion does not have the deficiencies of the Papapetrou solution [42ll43] in the 
original theory, and it is almost certainly indistinguishable from the Reissner- 
Nordstrom solution experimentally. Also, when this solution is expressed in 
Newman- Penrose tetrad form, it can be shown to be of Petrov Tvpe-D[24]. 
And of course the solution reduces to the Schwarzschild solution for Q — Q. 
And from ()115mi8| we see that the solution goes to the Reissner-Nordstrom 
solution exactly in the limit as vlf,— >cxd. 

The only significant difference between our electric monopole solution and 
the Reissner-Nordstrom solution occurs on the Planck scale. From (|108lll2p . 
the surface area of the solution is 1441. 



/ surface \ 
\ area / 



= / de 



2Q2 

Atr^ 



(121) 



The origin of the solution is where the surface area vanishes, so in our coordi- 
nates the origin is not at r = but rather at 



ro = yQ(2/^)i/4. 



(122) 
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From (|35|12p we have r^^lp ^ 10~'^^cm for an elementary charge, and tq <C 
2M for any reaUstic astrophysical black hole. For Q/M < 1 the behavior at 
the origin is hidden behind an event horizon nearly identical to that of the 
Reissner-Nordstrom solution. For Q/M > 1 where there is no event horizon, 
the behavior at the origin differs markedly from the simple naked singularity 
of the Reissner-Nordstrom solution. For the Reissner-Nordstrom solution all of 
the relevant fields have singularities at the origin, with goo^Q^ /r^ , Aq — Q/t, 
Foi=Q/i"^, -Roo~2(5'*/r^ and Rxi^2/r'^. For our solution the metric has a less 
severe singularity at the origin, with ^ — ^frj \Jt — tq. Also, the fields iV^,y, 
TV^''^ v^, A,, V^rt", ^U^, V^g"^, V^g.t., and the functions "a" 
and V all have finite nonzero values and derivatives at the origin, because it 
can be shown that f (ro) = ^/2 [^(l/4)]V6^/^-2/3 = f. 08137. The fields F^^, 
r"j^ and y/~gTZ^^ are also finite and nonzero at the origin, so if we use the 
tensor density form of the field equations (|41l45p . there is no ambiguity as to 
whether the field equations are satisfied at this location. 



9 An exact electromagnetic plane-wave solution 

Here we give an exact electromagnetic plane-wave solution for this theory 
which is identical to the electromagnetic plane-wave solution in Einstein- 
Maxwell theory, usually called the Baldwin-Jeffery solution[4F,46,47,48 . We 
will not do a full derivation, but a MAPLE program^^ which checks the solu- 
tion is available. We present the solution in the form of a pp-wave solution [17]. 
and a gravitational wave component is included for generality. The solution is 
expressed in terms of null coordinates x, y, u = (t — z)/\/2, v = {t + z)/-\/2. 




(123) 



(124) 



where 



fc^ = (0,0,-1,0), A^ = (0,0,AO), A = -V2{x,U+yfy), (125) 

H = 2H + A^ (126) 

= 2ih+x^ + h^xy - h+y^) + 2(/J + fl){x^ + y^), (127) 

H = h+x^ + hy,xy- h+y^ + {yf, - xfyf. (128) 



and the nonzero connections are 



33 - 2 9x ' 33 - 2 ay ' ^ 33 - 2 ay A, 



1 dH 2 d{fl + /2) 



du 
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IdH _ 2i dh ^,_ldH 2i dh 

23 2 9?; aw ' 32 2 dy ^ du' 

Here h^{u), (u) characterize the gravitational wave component, fx{u), fy{u) 
characterize the electromagnetic wave component, and all of these are arbitrary 
functions of the coordinate u = {t — z)/^/2. 

The solution above has been discussed extensively in the literature on 
Einstein-Maxwell theorv [151l461l471l45] so we will not interpret it further. It 
is the same solution which forms the incoming waves for the Bell-Szekeres 
colliding plane- wave solution [48], although the full BcU-Szekeres solution does 
not satisfy our theory because the electromagnetic field is not null after the 
collision. 



10 Conclusions 

The Einstein-Schrodinger theory is modified to include a cosmological constant 
Az which multiplies the symmetric metric. This is assumed to be nearly can- 
celled by Schrodinger's "bare" cosmological constant Ai, which multiplies the 
nonsymmetric fundamental tensor, such that the total cosmological constant 
A = Ab + Az matches measurement. The resulting theory closely approximates 
Einstein-Maxwell theory for jyl^j ^ 1/ (Planck length)'^, and it becomes exactly 
Einstein-Maxwell theory in the limit as >oo. 
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A Solution for iV,y^ in terms of g^^ and f^fj. 

Here we invert the definitions H4I22|I of and fun to obtain H33I34| I. the approximation of 
N^p, in terms of g^/^ and fvfi- First let us define the notation 

j^M=J^M^j^-l/2^ (130) 

We assume that If" fj,\^l for all components of the unitless field /"p, and find a solution 
in the form of a power series expansion in /"^p. Lowering an index on I I23I I gives 

(v^/V^)A^^^c = 5j:-r<.. (131) 

Let us consider the tensor f^^a = f^'^gva- Because g^a is symmetric and is anti- 

symmetric, it is clear that a =0. Also because fi^a-f" /j. is symmetric it is clear that 
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f" af^ iif^v = 0. In matrix language therefore tr(/)=0, tr(/^) = 0, and in fact tr{f^)=0 
for any odd p. Using the well known formula det(e^) = exp{tr{M)) and the power series 
ln{l-x) = -X- x^ 12 - a;^/3 - a;''/4 ... wo then getfio]. 



ln{detil-f)) = tr{ln{I-f)) = -/"..r p/2 + {/■*) . 



(132) 



Here the notation (/*) refers to terms like f'^af°'crf''pf''T- Taking ln{det()) on both sides 
of jmll using l(T32l l and the identities det{sM) = s'^detiM), det{M-^) = l/det(M) gives 



Taking on both sides of l|133ll and using = 1 + x + x^ /2 . . . gives 



1 



2(n-2) 



/""Up + (/*) ■ 



(133) 



(134) 



Using the power series (l — x) ^ = 1 + x + a;^ + a;^ . . or multiplying II131I I term by term, 
we can calculate the inverse of l|131|l to get |49) 



{V-<?/vMv)7v% = 5r, + /% + r,r„ + rpF^i^. + in ■ ■ ■ 



(135) 
(136) 



Here the notation (/*) refers to terms like fiyaf'^a-f'^pf^p- Since fu^f'^p is symmetric and 
fvpf''iyf'^p is antisymmetric, we obtain from II136I134I1301 the final result II33I34I I. 
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